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1.0 


INTRODUCTION 


^  This  report  summarizes  recent^ progress^ by  Science  Applications,  Inc., 
(SAI)  in  the  Investigation  of  low-density  air  chemistry  and  tensor  conductivity 
(Hall -current  effects)  and  its  Impact  on  high-current  beam  stability. 

a 

In  Section  ZxQrt  the  results  of  modeling  high  current  beams  in  low- 
density  air  are  discussed.  At  low-densities,  the  standard  simplifying  assumptions 
usually  employed  in  conductivity  modeling  no  longer  apply:  scalar  conductivity, 
ohm's  law;  local -Instantaneous  energy  deposition;  Maxwellian  distributed  plasma 
electrons;  no  delta  rays;  and  no  Inertial  effects.  These  assumptions  are  not 
made  in  the  present  model.  It  is  concluded  that  Hall  currents  do  play  a  signi¬ 
ficant  role  at  low  enough  densities  and  that  the  redistribution  of  plasma  cur¬ 
rent  can  result  In  a  significant  but  sudden  increase  in  the  magnetic  pinch  below 
a  "critical  ^ir  density. 

y  ‘ij 

Section  3.8-we  discuss-the  viability  of  using  the  AFWL/SNLA  HERMES 
facility  for  addressing  Important  Issues  of  high-current  hose  stability.  It  is 
concluded  that  with  proper  pulse  conditioning,  the  HERMES  pulse  should  be  able 
to  test  high-current  hose  growth  saturation  and  calibrate  existing  hose  codes. 


2.0 


HIQH  CURRENT  BEAMS  IN  LOW  DENSITY  AIR 


2.1  Introduction 


We  report  here  the  first  results  of  applying  a  phenomenological  low- 
density  chemistry  code  to  high  current  beams.  A  brief  account  of  the  model  has 
been  given  previously  (Refs.  1  and  2).  However,  substantial  modifications  were 
necessary  for  application  to  high  current  beams.  The  model  is  not  considered 
complete,  and  there. are  no  independent  calculations  in  the  same  parameter  regime 
to  which  It  can  be  compared;  thus  the  quantitative  results  are  tentative. 

Conclusions  based  on  the  calculations  are: 


1)  Below  a  certain  model -dependent  density  ^  .01  normal,  the 
electric  field  drives  a  bulk  runaway  which  changes  the  dis¬ 
tribution  of  plasma  currents.  The  net  result  is  a  sudden 
significant  Increase  In  the  pinch  force  as  the  density  is 
decreased  below  the  critical  value. 

2)  If  Hall  currents  are  turned  off,  the  enhanced  pinch  may  be 
reduced  by  as  much  as  a  factor  of  two  at  one  Bennett  radius. 

3)  Significant  amounts  of  plasma  current  are  driven  by  the 
gradient  in  electron  pressure. 

4)  Significant  amounts  of  plasma  current  are  carried  by  the 
high-energy  non-Maxwell  Ian  part  of  the  electron  distribution. 


This  section  Is  divided  Into  three  main  topics:  (1)  a  discussion 
of  the  physics  requirements  for  low-density  calculations,  (2)  an  outline  of 
the  present  status  of  the  model,  and  (3)  a  discussion  of  the  computational 
results  for  high  current  beans. 


2.2  Beam-Driven  Chemistry  In  the  Low-Density  Regime 

Simple  order-of-magnl tude  argments  show  that  at  background  densities 
sufficiently  lower  than  normal  atmospheric  density*  several  key  assumptions 
built  Into  standard  beam  chemistry  codes  are  not  justified.  Among  these  In¬ 
applicable  assumptions  are: 

1)  The  beam-initiated  cascade  results  In  local.  Instantaneous 
production  of  electrons  and  Ions. 

2)  Currents  associated  with  the  cascade  Itself  can  be  neglected. 

3)  Plasma  currents  can  be  calculated  from  Ohm's  law. 

4)  The  electron  and  Ion  densities  are  always  almost  equal,  so 
transport  effects  can  be  Ignored  In  the  chemistry  calculations. 

5)  The  electron  distribution  remains  close  to  Maxwellian,  so 
that  the  high  energy  parts  of  the  distribution  are  no  more 
Important  than  usual  In  determining  currents.  Ionization 
rates,  etc. 

In  addition,  many  beam  chemistry  codes  use  electromagnetic  algorithms 
which  Ignore  Hall  current  effects.  This  assumption  can  break  down  In  two  ways 
at  low  density:  (a)  the  momentum  transfer  frequency  goes  down  with  the  density, 
so  that  It  may  not  exceed  the  Larmor  frequency  by  a  wide  margin  as  In  full- 
density  air,  and  (b)  the  highly-overpopulated  high  energy  tall  of  the  electron 
distribution  at  low  density  may  have  a  momentum  transfer  frequency  considerably 
lower  than  that  of  the  bulk  of  electrons,  and  thus  make  a  larger  contribution 
to  the  current. 

In  the  past  (Ref.  3)  we  have  Investigated  the  effects  of  Hall  currents 
In  full  and  reduced  density  air,  but  have  not  been  able  to  address  problem  (b) 
above.  In  addition,  we  have  made  assumptions  1  -  5  In  our  chemistry  codes. 

At  present  we  are  developing  for  DARPA  a  multi -energy-group  model  which  abandons 
assumptions  1-5  and  addresses  the  non-Maxwell Ian  aspects  of  Hall  current  cal¬ 
culations.  Recent  calculations  by  Yu  (Ref.  4)  have  confirmed  the  Importance  of 
Hall  current  effects  for  ATA-Ilke  beams  when  the  electron  energy  distribution 
Is  treated  In  detail.  The  present  work  verifies  their  Importance  for  high  current 
beams. 
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: 


Non-Local,  Non-Ohmic  Conductivity  Model 


The  ultimate  goal  of  the  development  program  Is  to  produce  a  relatively 
simple  BMCOND-Ilke  model  which  can  be  Incorporated  Into  propagation  codes.  This 
requirement  eliminates  the  possibility  of  doing  a  "first  principles"  calculation. 
The  approach  taken  Is  to  develop  a  phenomenological  model  which  represents  the 
most  Important  physical  processes  in  a  simple  way,  and  to  adjust  the  "free" 
parameters  associated  with  various  simplifying  assumptions  to  obtain  agreement 
with  more  fundamentally-based  calculations. 


The  model  described  below  considers  the  electron  distribution  broken 
Into  three  energy  ranges.  The  lowest  energy  group  represents  basically  the  usual 
bulk  of  approximately-Maxwelllan  electrons.  The  next  higher  energy  group 
represents  those  that  In  the  presence  of  an  electric  field,  are  In  the  runaway 
regime  and  thus  behave  very  differently  from  the  bulk  of  the  electrons.  The 
third  group  represents  the  relativistic  particles  produced  directly  by  the  beam. 


The  organization  of  the  model  Is  summarized  in  Fig.  1,  In  which  the 
sources  and  sinks  of  particles  for  each  group  are  shown  schematically.  The 
three  energy  ranges  will  be  referred  to  by  the  names  low-group,  high-group, 
and  6 -group  In  order  of  increasing  energy.  The  beam  particles  themselves 
constitute  a  fourth  group  which  Is  treated  In  the  usual  way  (e.g.,  no  energy 
straggling). 


The  beam  particles  colllslonally  produce  secondaries  directly  In  each 
of  the  three  groups,  according  to  a  Moller  distribution  (modified  at  low  energy). 
The  maximum  5-ray  energy  is  one-half  the  beam  energy,  since  the  higher  energy 
particle  emerging  from  a  primary  Interaction  remains  associated  with  the  beam. 
Collisions  by  beam  particles  are  the  only  source  of  electrons  In  the  6 -group. 

A  more  detailed  discussion  of  the  5-ray  model  Is  given  below  In  Section  2.3.1. 
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2.1.  3-Group  Model 


The  high-energy  group  extends  downward  from  1  MeV  to  a  variable-energy  boundary 
which  at  present  goes  no  higher  than  10  KeV.  It  is  fed  mainly  by  primary  beam 
collisions,  and  by  the  colllslonal  degradation  of  the  5-group  through  Its  1  MeV 
lower  boundary.  Under  some  conditions,  electrons  are  also  transferred  up  from 
the  low  group. 

Both  the  high  and  low  energy  electron  groups  are  represented  by  fluid 
equations  for  particle  density,  momentum  and  energy.  The  equations  of  motion 
provide  the  means  for  abandoning  the  ohmic  representation  of  plasma  current 
flow.  The  fluid  equations  are  solved  on  an  Eulerian  grid  which  is  the  same 
grid  used  for  the  electromagnetic  field  calculation.  It  is  assumed  that  every¬ 
thing  is  a  function  of  the  "retarded"  time  variable  £  =  t  -  z/c  only,  so  that 
the  only  independent  variables  are  r,  the  radial  position,  and  £»  the  distance 
back  from  the  pulse  head  in  seconds.  The  electromagnetic  fields  are  obtained 
using  the  algorithm  described  in  Ref.  3. 

In  principle  it  is  necessary  to  solve  the  fluid  equations  and  5 -group 
equations  simultaneously  with  the  electromagnetic  field  equations.  Since  the 
Maxwell  equations  themselves  are  explicitly  non-linear  when  the  Hall  current 
terms  are  Included,  and  the  fluid  equations  are  also  non-linear,  an  iterative 
procedure  is  required.  In  order  to  maintain  maximum  modularity  under  these 
circumstances,  the  procedure  adopted  involved  solving  first  a  linearized  set 
of  fluid  equations,  followed  by  the  field  equations,  and  then  iterating  to 
convergence. 

In  addition,  the  algorithm  has  been  designed  to  allow  differential 
comparisons  between  the  three-group  model,  and  a  model  with  the  high  energy 
and  5-groups  eliminated  and  with  the  standard  local-instantaneous  approximations 
made.  However,  this  latter  "one-group  model"  still  differs  from  standard 
chemistry  codes  in  that  it  solves  an  equation  of  motion  for  the  one  group 
rather  than  using  Ohm's  law.  It  is  also  possible  to  turn  off  the  pressure 
terms  in  the  equations  of  motion,  and  the  Hall  current  effects. 


2.3.1 


Delta-Ray  Group 


In  the  present  discussion,  6-rays  will  refer  to  the  part  of  the  beam- 
initiated  cascade  which  has  energy  greater  than  1  MeV.  These  particles  are 
represented  by  a  model  essentially  identical  to  that  developed  by  Johnston 
(Ref.  5).  As  applied  here,  the  model  has  two  main  deficiencies:  (1)  It  does 
not  provide  detailed  information  on  the  radial  distribution  and  (2)  it  does  not  - 
include  the  effects  of* electric  or  magnetic  fields  on  the  5-particles.  With 
the  assumption  that  everything  depends  on  £  =  t  -  z/c,  the  model  gives  an 
*itegral  expression  for  the  total  6-ray  current  as  a  function  of  distance  £ 
back  from  the  beam  head.  Particle  energy  distribution  information  is  available 
but  not  used  at  present. 

It  Is  assumed  that  most  of  the  particles  produced  as  secondaries  by 
the  6-rays  have  very  low  energy.  This  Is  consistent  with  the  continuous  slowing 
down  approximation  which  forms  the  basis  of  the  model.  In  the  present  model, 
as  the  6-rays  lose  energy,  they  produce  low  energy  electrons  at  a  rate  of  one 
particle  per  33.73  eV  of  energy  lost;  these  particles  immediately  enter  the  low- 
energy  plasma  electron  group  of  the  three-group  model.  The  6-rays  themselves 
degrade  in  energy  due  to  collisions.  They  reach  the  cutoff  of  1  MeV  at  a  rate 
given  by  the  Johnston  model,  and  at  this  point  they  are  added  to  the  high  energy 
electron  group. 

The  radial  profile  as  a  function  of  distance  back  from  the  beam  head 
is  estimated  in  a  crude  way,  based  on  two  pieces  of  information:  (1)  the  5's 
are  produced  with  the  beam  profile  and  (2)  they  evolve  from  that  profile  to 
form  a  halo  with  radial  dimension  estimated  by  Johnston  (Ref.  5)  to  be  ^  4 
Bennett  radii  about  a  Bennett- profile  beam  (for  particles  above  1  MeV).  It 
is  assumed  that  at  a  given  £,  the  profile  towards  which  the  5*s  evolve  is  a 
weighted  average  of  the  beam  profile  and  a  specified  halo  profile.  The  relative 
weights  are  determined  from  the  rate  of  local  production  as  compared  to  the 
total  Instantaneous  number  density,  and  the  rate  of  evolution  towards  the  halo 
profile.  This  assumption  is  used  in  the  following  interpolation  formula  for 
the  radial  profile  f(r,  C+b£)  in  terms  of  f(r,£): 


f(r,  C+AC)  »  (f(r,C)  +  A$F)/(1  +  vdAC). 

The  function  F  is  given  by 

F  *  (§  h  +  g  v), 

where  S  is  the  local  rate  of  5-ray  production  by  the  beam,  N  is  the  total  number 
of  5-rays  present,  and  h  is  the  beam  profile,  on  which  the  5-rays  are  produced. 
The  profile  g  is  the  assumed  halo  profile,  and  v  is  a  characteristic  rate  at 
which  the  5's  evolve  towards  the  profile  g.  The  rate  is  set  equal*  to  S/N  +  v. 
With  this  choice,  f(r,£+AC)  is  normalized  to  unity  if  h,  g,  and  f(r,£)  are. 

The  rate  v  is  taken  to  be  an  average  radial  velocity  of  the  5's  divided  by  the 
assumed  halo  characteristic  radius. 

2.3.2  The  High-Energy  Group 

The  schematic  forms  of  the  fluid  equations  for  the  high-group  are: 

||  +  V-(Nv)  +  VjN  =  SB  +  S5  +  SL  (2.1) 

If  +  7*<3?)  +  i  t  v2J  -  ^  \  &  (f  ♦  i  XB)  (2.2) 

~  (|  Nmv2  +  |  P)  +  V«[v(|  Nmv2  +  |  P)]  +  V*(vP) 

+  Vj(  \  Nmv2  +  |  P)  «  L  +  NeE*v  (2.3) 

The  terms  Sg,  Kg,  etc.  are  the  sources  due  to  the  beam,  the  5-group,  and  upward 
transfer  from  the  low-group,  respectively.  It  is  assumed  that  all  collisional 
ionizations  by  high-group  electrons  produce  a  low  energy  secondary  which  goes 
directly  into  the  low  group.  The  energy  equation  (2.3)  has  been  written  with 
the  total  energy  split  into  a  drift  part,  and  a  thermal  part  represented  by  a 
pressure  P.  The  term  L  Includes  energy  deposited  by  the  beam,  energy  brought  In 
by  transfer  of  particles  from  the  low-group  and  5-group,  and  energy  losses  due 
to  ionization.  Energy  loss  due  to  loss  of  particles  from  the  high  group  is 


represented  by  the  term  on  the  left  side  of  the  equation.  It  has  been  assumed 
that  the  particles  lost  take  with  them  the  average  energy  of  the  high  group.  The 
momentum  loss  rate  v2  Includes  the  effect  of  particle  losses,  momentum  loss  due 
to  Ionizing  collisions,  and  momentum  loss  due  to  elastic  collisions. 

The  fluid  equations  are  solved  as  a  finite  difference  systemonthe  Eulerlan 
electromagnetic  field  grid.  The  outer  boundary  condition  allows  an  outflow  of 
particles.  No  explicit  assumption  is  made  concerning  the  distribution  function 
of  the  high  group  electrons  except  at  5  *  0+,  when  it  is  dominated  by  the  Moller 
distribution  from  primary  beam  particle  Interactions.  The  average  total  energy 
and  z-momentum  are  then  used  to  give  an  initial  "temperature"  T,  defined  by  P  *  NT, 
where  P  is  the  pressure  from  Eq.  (2.3).  At  all  subsequent  times,  the  temperature, 
particle  density,  and  drift  velocity  are  obtained  from  the  fluid  equations. 

For  simplicity,  ionization  and  momentum  transfer  rates  are  evaluated 
at  the  mean  particle  energy.  An  analytic  formula  given  by  Briggs  and  Yu  (Ref.  6) 
is  used  for  the  ionization  cross  section.  The  momentum  transfer  frequency  vm  Is 
given  by  the  approximate  formula 


=  (1.08  x  10"5  N  &*)/(178.89  +&/Z) 

g 

+  (1.46  x  10“6  N+  log  A)/T^iv 


(2.4) 


where  £  is  the  total  electron  energy  in  volts.  The  first  term  represents  the 

effect  of  collisions  with  neutral  particles,  and  the  second  represents  Coulomb 

2 

collisions.  The  "equivalent"  temperature  is  defined  as  T  *  T  +  1/3  mv  , 

log  A  is  the  usual  Coulomb  logarithm,  Na  is  the  total  density  of  neutral  atoms 

+  * 

and  molecules,  and  N  is  the  total  density  of  positive  ions. 

The  boundary  between  the  high  and  low  energy  groups  is  set  by  finding 
the  solution  of  the  equation 


e  E  *  m  v 
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corresponding  to  the  high  energy  side  of  the  peak  In  The  average  velocity 
v  Is  related  to  the  total  energy  £  appearing  In  Eq.  (2.4)  byS  ■  %  m  v2. 

However,  the  boundary  energy  Is  not  allowed  to  exceed  10  KeV  even  when  the 
total  electric  field  E  Is  very  small. 

Another  special  case  occurs  when  the  field  Is  so  strong  that  there  Is 
no  physical  solution  to  (2.5);  this  corresponds  to  a  bulk  runaway.  In  which  the 
field  Is  capable  of  accelerating  the  electrons  through  the  peak  of  the  momentum 
transfer  cross  section.  In  this  case,  the  choice  of  the  lower  boundary  of  the 
high  group  is  somewhat  arbitrary.  The  lower  limit  to  be  used  for  the  boundary 
during  a  strong  runaway  condition  can  be  specified  as  Input  data.  Usually  It 
Is  taken  to  be  ^  a  few  volts,  so  that  essentially  all  the  electrons  produced 
directly  by  the  beam  go. Into  the  low  end  of  the  high  energy  group.  The  transfer 
of  electrons  from  the  low  to  the  high  group  Is  discussed  below. 

2.3.3  The  Low-Energy  Group 

The  equations  for  the  low  energy  group  are  similar  In  form  to  those 
given  above  for  the  high  group  with  the  exception  that  they  contain  loss  terms 
due  to  recombination  and  to  transfers  to  the  high  energy  group,  and  input  terms 
due  to  transfers  from  the  high  group  and  to  collisional  ionization  by  high 
group  electrons.  The  main  difference  in  the  treatment  of  the  low  group  Is  In 
the  chemistry  detail.  At  present,  the  abundance  of  Ng,  02,  N,  N(  D),  0,  and  a 
composite  representative  of  the  triplet  states  of  N2  are  followed  explicitly  by 
differential  equations.  The  system  of  equations  used  Is  essentially  the  same 
as  In  the  chemistry  code  BMCOND  (Ref.  8).  This  degree  of  complexity  was  found 
necessary  to  obtain  reasonable  agreement  with  the  comprehensive  code  HICHEM  at 
low  air  densities  for  high  current  calculations. 


Ionization  rates  are  given  by  tables  Identical  to  those  in  the  HICHEM 
code  up  to  500  volts  "equivalent"  temperature  (TeqUjv  =  Te  +  1/3  m  v  ).  The 
collisional  excitation  rates  from  the  8MC0ND  model  are  used.  Ionization  and 
excitation  rates  due  to  the  beam  are  computed  as  In  the  HICHEM  code,  but  with 
appropriate  adjustment  for  the  differences  between  time-delayed  deposition  model 
used  here  and  the  Instantaneous  deposition  model  In  HICHEM  and  BMCOND. 

It  Is  very  important  to  distinguish  between  the  positive  ion  density 
and  the  electron  density,  since  electron  transport  effects  can  be  very  large; 
thus  a  separate  differential  equation  for  the  total  rate  of  positive  ion  produc¬ 
tion  Is  integrated.  In  addition  It  was  found  necessary  to  calculate  a  vibrational 
temperature  because  of  significant  sensitivity  of  computational  results  to  the 
dissociative  recombination  rate;  the  HICHEM  treatment  was  adopted  for  this 
calculation.  The  momentum  transfer  frequency  used  In  BMCOND  Is  used  for  low 
energies,  with  a  smooth  transition  to  the  analytic  formula  (Eq.  2.4)  used  for 
the  high  group.  The  Coulomb  term  of  Eq.  (2.4)  is  used  throughout. 

Reassignment  of  electrons  from  the  low  group  to  the  high  group  Is  a 
relatively  arbitrary  procedure.  The  goal  Is  to  remove  electrons  from  a  presumed 
high-energy  non-Maxwell Ian  tall  of  the  low  group  and  put  them  Into  the  high 
group,  which  hopefully  represents  better  their  contribution  to  currents,  ioniza¬ 
tion  rates,  etc.  than  does  the  assumed  Maxwellian  bulk  of  the  low  group  electrons. 
However,  there  is  no  simple  model  of  the  super-thermal  tail  as  generated  by 
strong  electric  fields  which  vary  rapidly  in  both  time  and  space.  Several 
different  ad  hoc  procedures  have  been  tried,  but  none  is  especially  defensible 
In  detail. 

Qualitatively,  what  happens  in  a  weak  field  (or  high  ambient  density 
gas)  Is  not  very  sensitive  to  the  details  of  the  transfer  rate.  Under  strong 
runaway  conditions,  the  entire  low  group  is  accelerated  to  high  drift  energy 
and  heated  to  temperatures  ^  kilovolts  on  a  very  short  timescale.  During  this 
time,  the  distinction  between  the  two  groups  is  not  very  meaningful,  and  again 
the  behavior  of  the  bulk  plasma  is  not  too  sensitive  to  the  details  of  the  transfer 


rate.  However,  after  the  fields  which  precipitate  a  strong  runaway  die  out, 
the  distribution  of  electrons  between  the  two  groups  can  be  important;  this 
does  depend  strongly  on  the  transfer  rate.  Similarly,  If  conditions  for  a 
strong  runaway  are  just  barely  achieved,  for  only  a  short  time,  there  may  be 
some  sensitivity  to  the  transfer  rate.  These  problems  are  continually  being 
studied,  both  In  the  context  of  the  model  described  here,  and  through  comparison 
with  other,  more  comprehensive  models  being  developed. 

The  transfer  rate  presently  used  Is  calculated  as  follows.  A  velocity 
v  isdeflned  by  v„  *  +  v  +  .5  Vt,  where  the  +  sign  Is  taken  If  the  radial  drift 
velocity  Is  parallel  to  Er,  and  a  negative  sign  Is  taken  If  It  is  opposed.  The 
thermal  velocity  term  is  added  to  account  for  the  fact  that  at  high  enough 
temperature,  a  substantial  number  of  electrons  may  be  able  to  run  away  even 
If  the  bulk  drift  Is  not  large  (or  even  parallel  to  the  Er  field).  The  transfer 
rate  of  electrons  Is  then  given  by 

v  -v 

SL  =  min(.5,  .lNh/Ne)  min(l.,  exp{-\±))/HZ ,  (2.6) 

in  which  v  Is  the  velocity  corresponding  to  the  energy  boundary  between  the 
groups,  Nh  is  the  local  density  of  high-group  electrons,  Ne  Is  the  local  density 
of  low-group  electrons,  and  is  the  proposed  tlmestep.  The  purpose  of  the  first 
minimum  function  Is  to  assure  that  the  high-group  density  is  not  changed  by  more 
than  10%  during  the  step.  In  actual  computations,  the  transfer  rate  is  usually 
very  low,  or  else  as  high  as  permitted  by  the  first  factor  of  (2.6),  even  for 
exceedingly  small  tlmesteps. 

The  momentum  transferred  with  the  particle  is  taken  to  be  parallel 
to  Its  total  drift  velocity,  with  magnitude  given  by  the  larger  of  the  velocities 
v  and  v  .  The  energy  transferred  is  the  average  energy  per  particle  of  the 

v  C 

low  group,  plus  an  additional  amount  arranged  to  come  (by  suitable  terms  In 
the  equations  of  motion)  entirely  from  the  drift  energy  of  the  low  group,  to 
make  up  the  total  drift  energy  of  the  electron  injected  into  the  high  group. 


12 


The  net  result  Is  to  decrease  the  drift  energy  of  the  low  group,  but  leave  Its 
temperature  unchanged.  The  effect  of  the  energy  and  momentum  transfers  on  the 
high  group  depends  on  its  energy  and  momentum  at  the  time  of  the  transfer. 

Other  procedures  for  transferring  electrons  to  the  high  group  are  being  evaluated. 

2.4  Calibration  in  the  High  Current  Regime 

At  present  there  are  no  comprehensive,  first-principles  calculations 
to  provide  detailed  quidance  in  the  development  of  the  three-group  model  for 
high  current  calculations.  Yu  (Ref.  4)  has  presented  Boltzmann-code  calcula¬ 
tions  for  ATA-like  parameters,  but  even  these  calculations  incorporate  assump¬ 
tions  which  are  not  always  justified  in  the  context  of  high  power  beam  plasmas. 
The  purpose  of  the  calibration  runs  discussed  below  is  to  check  that  the 
relatively  simple  model  described  in  Section  2.3  agrees  reasonably  well  with 
detailed  chemistry  codes  in  regimes  where  the  assumptions  of  those  codes  (see 
Section  2.2)  are  not  thought  to  be  seriously  in  error.  Results  of  calibration 
at  10  kA  have  been  presented  elsewhere  (Ref.  1).  A  similar  comparison  with  the 
HICHEM  (Ref.  7)  code  at  100  kA  and  gas  density  0.1  normal  is  given  below.  The 
rise  time  is  5  ns  and  the  Bennett  radius  is  0.5  cm. 

The  comparison  of  electron  densities  on  axis  and  at  one  Bennett  radius 
is  shown  in  Figure  2.2.  The  agreement  is  very  satisfactory.  The  effect  on  the 
electron  density  due  to  the  time  delay  for  the  three-group  model  is  not  large 
after  2  or  3  ns.  The  electron  density  on  axis  is  mainly  determined  by  the 
close  balance  between  collisional  ionization  and  dissociative  recombination 
after  a  few  ns.  The  recombination  rate  itself  decreases  significantly  as  the 
molecules  (and  molecular  Ions)  are  depleted,  and  as  the  vibrational  temperature 
increases.  After  about  6  ns,  the  continued  increase  in  Ng  is  determined  mainly 
by  the  decrease  in  the  recombination  rate,  with  a  large  part  attributable  to 
the  vibrational  temperature  dependence. 


A  similar  comparison  of  electron  temperatures  and  total  currents  Is 
shown  in  Figure  2.3.  The  temperature  agreement  is  quite  good  from  1  ns  on. 

At  0.1  ns,  HICHEM  gives  Tg  'v  18  volts,  compared  to  the  12  or  13  volts  given 
by  the  simple  models.  At  such  early  times  the  ohmic  approximation  used  In 
HICHEM  is  not  very  good,  so  the  ohmic  heating  rate  is  likely  to  be  Incorrect; 
thus  It  is  not  clear  which  values  are  closer  to  the  truth. 

As  in  the  10  kA  comparison  (Ref.  1),  the  agreement  between  HICHEM 
and  the  simple  models  is  not  as  good  for  the  net  current  and  effective  current. 
This  difference  is  mainly  due  to  the  difficulty  in  matching  the  momentum 
transfer  frequency  calculated  in  HICHEM  by  a  very  simple  formula.  However, 
the  agreement  in  currents  is  still  acceptable,  particularly  since  the  most 
Interesting  sensitivities  described  below  (Section  2.5)  develop  before  1  ns. 

Because  the  results  of  the  HICHEM  code  are  suspect  below  p/p0  %  0.1, 
especially  for  £  <  f ew  ns,  it  cannot  be  used  for  calibration  comparisons  at 
lower  densities.  However,  many  of  the  results  presented  below  are  differential 
comparisons  or  sensitivity  studies,  and  thus  can  provide  useful  Information  in 
spite  of  uncertainties  in  the  quantitative  results. 

2.5  Computational  Results 

2.5.1  Introduction 

The  beam  parameters  for  the  calculations  described  below  are: 

Current  *  100  kA 

Rise  time  *  5  ns 

Bennett  radius  »  0.5  cm 
Energy  »  10  MeV 
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The  discussion  below  Is  divided  Into  two  main  sections.  The  sen¬ 
sitivities  to  gas  density.  Hall  currents,  and  pressure  terms  are  described 
briefly  In  Section  2.5.2,  and  the  remaining  sections  discuss  the  calculations 
In  more  (detail. 

The  main  conclusions  have  been  stated  earlier  (Section  2.1).  In 
brief.  It  has  been  found  that  the  pinch  force  is  strongly  dependent  on  the  gas 
density  below  some  threshold  value,  and  significantly  dependent  on  the  presence 
of  Hall  currents  and  pressure  terms.  The  plasma  current  and  electron  density 
radial  distributions  change  dramatically  over  a  very  small  range  of  density 
ratio  near  p/pQ  *  .01,  resulting  In  an  Increase  of  pinch  force  by  a  factor 
of  3  or  more  for  a  density  change  of  only  ^  20%. 

2.5.2  Sensitivity  Results 

2.5.2. 1  Density  Sensitivity  of  Pinch  Force 

The  dependence  of  I@ff(r)  on  p/pQ  Is  shown  In  Figure  2.4  for  the 
three-group  model,  and  In  Figure  2.5  for  the  one-group  model.  The  sharp  onset 
of  current  enhancement  begins  In  the  range  p/pQ  ^  .0085  -  .01  for  the  three- 
group  model,  and  between  .005  -  .008  for  the  one-group  model.  In  both  cases, 
the  current  enhancement  Is  a  factor  ^  3.5  or  greater,  compared  with  a  base 
level  at  p/pQ  ^  .01.  These  large  effects  set  In  when  the  plasma  electrons  can 
be  sustained  In  a  state  of  bulk  runaway  for  several  tenths  of  a  nanosecond. 

The  sharp  density  threshold  results  from  the  sensitivity  of  the  charge 
neutralization  process  to  gas  density.  Several  effects  contribute.  The  beam 
production  of  positive  Ions  Is  directly  proportional  to  the  density;  the 
electrons  produced  by  the  beam  provide  "seeds"  for  the  avalanche  process  of 
plasma  electron  collisions.  The  peak  e-foldlng  rate  of  the  avalanche  is  also 
proportional  to  the  density.  These  two  processes  clearly  delay  the  charge 
neutralization  as  the  gas  density  is  decreased  and  allow  higher  radial  electric 
fields  to  develop  (assuming  the  rise  time  Is  not  «  1  ns).  In  order  for 
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neutralization  to  occur,  plasma  electrons  produced  by  various  processes  must 
move  out  of  the  spatial  region  occupied  by  the  beam,  leaving  the  positive  Ions 
to  cancel  the  charge  of  the  beam  particles.  Higher  fields  help  move  the 
electrons  out  quickly,  but  that  process  dilutes  very  considerably  the  avalanche 
Ionization  and  thus  slows  the  production  of  the  needed  positive  Ions  close  to 
the  beam  axis.  In  the  true  IFR  regime,  avalanche  is  relatively  unimportant,  and 
the  accumulation  rate  of  beam-generated  positive  Ions  determines  the  neutral¬ 
ization  time  and  the  peak  fields. 

The  difference  In  density  threshold  between  the  two  models  is  caused 
by  (1)  the  time-delay  In  Ionization  in  the  three-group  model  which  reduces  the 
effective  beam  Ionization  rate  by  more  than  a  factor  of  two;  and  (2)  the  effects 
of  the  non-Maxwell Ian  high-energy  group  on  the  ionization  rate  and  on  the 
movement  of  plasma  electrons  away  from  the  beam.  It  is  difficult  to  assess 
these  separately  because  beam  production  of  the  high-energy  and  5 -ray  groups 
(which  degrade  relatively  slowly  at  low  densities)  is  the  cause  of  the  time  • 
delay.  Since  qualitatively-similar  results  occur  whether  or  not  the  high- 
energy  group  Is  Included,  it  seems  that  the  time-delay  Is  probably  most  directly 
responsible. 

2. 5. 2. 2  Sensitivity  to  Hall  Currents  and  Electron  Pressure 

The  radial  pressure  gradient  in  the  equations  of  motion  for  the  high 
and  low  energy  electrons  can  drive  currents  both  radially  and  In  the  z-dlrection 
due  to  the  magnetic  part  of  the  Lorentz  force  (Hall  effect).  This  works  in 
conjunction  with  the  electrically  driven  currents.  The  result  of  deleting 
either  the  pressure  terms  or  all  Hall  effects  (by  zeroing  the  magnetic  force 
on  the  plasma  electrons)  Is  shown  In  Figure  2.6.  Clearly,  the  Hall  currents 
have  the  largest  effect,  but  the  pressure  terms  are  not  negligible  either. 

Note  that  these  comparisons  show  the  cumulative  effect  of  removing  the  terms  • 
for  the  entire  calculation,  and  not  simply  the  contributions  to  the  current 
at  the  time  shown.  A  similar  calculation  was  done  with  only  the  high  group 
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pressure  turned  off.  The  effect  at  £  ■  1  ns  was  considerably  smaller  than  shown 
In  the  figure.  Low-group  pressure  contributes  dominantly  during  the  bulk 
runaway  because  the  electrons  reach  temperatures  In  excess  of  a  kilovolt,  and 
their  number  density  Is  large.  However,  It  will  be  shown  In  later  discussion 
that  the  high-group  pressure-driven  Hall  current  Is  Important. 

2.S.3  Comparison  of  1-Group  and  3-Group  Calculations  at  p/pQ  »  .008 

The  purpose  of  the  following  discussion  Is  to  provide  a  more  detailed 
description  of  the  phenomena  which  lead  to  the  enhanced  pinch  force.  At 
p/pQ  *  .008,  the  3-group  model  shows  a  very  strong  effect,  whereas  the  1-group 
model  shows  very  little  because  Its  density  threshold  Is  somewhat  lower. 

2.5. 3.1  Effects  on  Plasma  Current  Distribution 

The  net  current  (including  displacement)  integrated  out  to  radius  r 
Is  given  by  Inet  *  .005  r  B(r),  where  I  Is  In  kA,  B  Is  In  gauss,  and  r  Is  in  cm. 
The  effective  current,  which  measures  the  pinch  force,  is  the  beam-profile- 
weighted  average  of  .01  r  (B(r)  -  Er(r)).  The  value  of  the  effective  current 
Integral  taken  to  radius  r,  and  the  net  current,  are  shown  in  Figures  2.7(a)  -  (h) 
at  various  distances  from  the  pulse  head.  Large  differences  are  apparent  by 
1  ns  and  persist  to  10  ns. 

At  0.1  ns,  the  net  currents  for  the  two  calculations  are  similar, 
but  the  effective  current  Is  substantially  weaker  for  the  3-group  case  because 
the  radial  electric  field  is  higher  (due  to  slower  charge  neutralization). 

At  1  ns  the  effective  current  In  the  3-group  model  is  'v  28%  of  the 
beam  current  at  that  time,  whereas  in  the  1-group  model  It  is  only  ^  7.5%. 

The  net  current  profiles  Imply  a  very  much  broader  plasma  current  distribution 
for  the  3-group  model,  although  net  currents  inside  5  cm  radius  (10  Bennett 
radii)  are  ^  .5  to  1  kA  In  both  cases  even  at  10  ns  Into  the  pulse  where  the 
beam  current  Is  almost  100  kA. 


A  comparison  of  the  plasma  current  (Jz)  radial  profiles  at  1  ns  Is 
shown  In  Figure  2.8.  The  breakdown  of  the  3-group  profile  into  contributions 
from  the  high  and  low-energy  electron  groups  is  discussed  below  in  Section  2. 5. 4. 4. 

2. 5. 3. 2  Electron  and  Ion  Density  Profiles 

Major  qualitative  differences  in  the  evolution  of  the  electron  density 
are  apparent  in  Figures  2.9(a)  -  (d).  The  1-group  calculation  shown  in 
Figure  2.9(a)  has  the  usual  peak  slightly  off  axis  due  to  the  Er-initiated 
avalanche,  and  is  not  unusual.  The  3-group  calculation  for  the  low  group 
(Figure  2.9(b))  is  much  more  interesting.  It  shows  a  large  peak  well  off 
axis,  and  much  higher  density  at  large  radius.  The  high-group  density 
(Figure  2.9(c))  has  features  which  invite  interpretation  as  propagating  sound 
waves.  The  ridge  which  appears  earliest  in  time  is  associated  with  the  rise 
of  the  Er  field,  and  the  second  begins  on  the  axis  near  the  peak  of  Ez.  There 
is  a  hint  of  the  first  ridge  in  the  low-group  density  also. 

The  total  positive  ion  density  is  shown  in  Figure  2.9(d).  Since 
the  immobile  ions  show  the  same  gross  features  as  the  low-group  electron 
density,  it  is  clear  that  the  large  off-axis  peak  must  be  interpreted  in  terms 
of  the  history  of  the  ionization  rates.  However,  the  weak  ridge  in  the  low- 
group  electron  density  does  not  appear  in  the  ion  density,  and  thus  may  be  a 
flow  feature.  The  number  of  particles  in  the  high  group  is  not  large  enough 
for  the  ridges  of  Figure  2.9(c)  to  show  up  in  the  ion  density  (Figure  2.9(d)). 

It  seems  likely  that  they  are  similar  in  nature  to  the  ridge  in  the  low-group 
electron  density.  Further  support  for  the  flow  explanation  is  provided  by 
Figures  2.10(a)  and  (b).  These  show  the  low-group  and  total  ion-densities 
for  the  1-group  calculation  at  p/pQ  =  .005,  in  which  there  is  a  large  current 
enhancement  compared  to  p/pQ  =  .008  (see  Figure  2.7).  Here  the  ridge 
structures  in  the  electron  density  are  more  prominent,  but  still  have  no  strong 
counterpart  in  the  ion  density. 


The  difference  In  positive  Ion  density  radial  profiles  at  1  ns  Is 
shown  clearly  In  Figure  2.11(a).  At  this  time  the  Ion  density  Is  much  greater 
than  the  beam  particle  density,  so  these  profiles  are  very  close  to  the  electron 
density  In  the  charge-neutralized  state.  The  effect  on  the  z-component  of 
plasma  current  was  shown  earlier  (Figure  2.8).  Figure  2.11(b)  shows  how  the 
relatively  small  difference  In  Ion  densities  on  axis  due  to  time  delay  before 
0.1  ns,  becomes  very  large  for  a  few  tenths  of  ns  before  coming  together  again. 


For  ease  of  detailed  comparison,  electron  density  radial  profiles 
for  the  1-group  model,  and  for  the  low  and  high-energy  components  of  the  3- 
group  model,  are  shown  in  Figures  2.12(a)  -  (c).  The  large  temporary  off-axis 
hump  In  the  low-group  electrons  Is  very  clear  at  0.6  ns  In  Figure  2.12(b). 

The  very  much  broader  density  profiles  produced  by  the  3-group  model  are  also 
very  obvious,  and  persist  out  to  10  ns  from  the  pulse  head. 


A  summary  of  electron  and  Ion  number  density  comparisons  on  axis  Is 
given  in  Table  2.1  below.  The  first  two  entries  of  the  third  column  show  the 
magnitude  of  the  time-delay  effect,  while  subsequent  entries  show  the  large 
differences  seen  in  Figure  2.11(b).  The  first  two  columns  give  an  indication 
of  the  importance  of  electron  transport.  At  0.3  ns  for  the  3-group  model,  only 
47%  of  the  plasma  electrons  ever  produced  on  axis  remain  there. 


TABLE  2.1 

ELECTRON  TRANSPORT  EFFECTS 


1-Group 


3-Group 


.69 

.68 

.37 

.014 

.027 

1.27 

.82 
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Figure  2.11(a).  Positive  Ion  Radial  Profile 
at  1  ns  (p/pQ  *  .008). 


Figure  2.11(b).  Positive  Ion  Density  on  Axis 


Fig.  2.12(a).  Low  Energy  Group  Electron  Fig.  2.12(b).  Low  Energy  Group  Electron 

Density  Profiles,  1-Group  Density  Profiles,  3-Group 

Model  (p/pQ  *  .008).  Model  (p/P()  »  .008). 


R  (cm) 


Fig.  2.12(c).  High  Energy  Electron  Density 
Profiles,  3-Group  Model 


2.5.3. 3  Electric  Field  Effects 

Both  the  Ez  spike  and  the  final  decay  of  Ef  occur  significantly  earlier 
In  time  for  the  1-group  model  (Figures  2.13(a)  -  (d)),  and  the  peak  field  strengths 
are  lower.  The  large  negative  oscillation  In  Ep  at  1  Bennett  radius  Is  especially 
prominent  In  the  3-group  calculation,  but  appears  to  damp  out  satisfactorily. 

The  negative  spike  In  Er  In  caused  by  an  overshoot  of  the  outward-moving  electrons 
as  charge  neutralization  Is  finally  achieved  near  the  axis.  The  significant 
overshoot  Is  consistent  with  the  fact  that  the  electron  plasma  frequency  Is  much 
greater  than  the  collision  frequency  at  a,  0.3  ns.  The  large  off-axis  hump  In 
Ion  density  at  0.6  ns  (Figure  2.12(b))  seems  to  be  associated  with  the  decelera¬ 
tion  of  the  low-group  electrons  by  the  reversed  Er,  and  by  their  final  cooling 
through  the  peak  of  the  Ionization  cross  section.  The  final  cooling  occurs 
considerably  earlier  off-axis,  as  shown  In  Figure  2.14(a).  The  temperature 
behavior  of  the  1-group  model  Is  also  shown  for  completeness  In  Figure  2.14(b). 

2.5.4  Hall  Current  Effects 

2.5.4. 1  Introduction 

Hall  currents  may  be  comparable  to  ordinary  currents  when  the  Larmor 
frequency  eB/mc  Is  comparable  to  or  greater  than  the  momentum  transfer  frequency. 
This  criterion  Is  a  strong  function  of  both  ambient  density  and  electron  energy. 

The  3-group  model  allows  the  possibility  of  accounting  for  the  energy  dependence 
in  a  strongly  non-Maxwell Ian  plasma.  The  momentum  transfer  frequency  at  1  eV 
In  N2  Is  about  the  same  as  that  at  200  eV,  and  decreases  with  energy  beyond  50 
or  60  eV.  Thus,  If  a  substantial  fraction  of  the  electrons  have  energy  in  excess 
of  200  eV,  their  low  momentum  transfer  frequency  may  result  In  a  significant  Hall 
current,  out  of  proportion  to  their  fractional  number  density. 
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2. 5.4.2  Detailed  Discussion 


It  Is  clear  from  Figure  2.6  discussed  In  Section  2.3  above  that  Hall 
currents  are  Important  to  the  pinch  force.  In  the  calculation  at  p/pQ  »  .008, 
even  the  low-energy  electron  group  achieves  high  velocity  for  a  short  time, 
and  contributes  substantial  Hall  current.  However,  this  Is  generally  masked 
by  the  Ez-driven  return  current.  In  the  far  wings,  however,  and  at  early 
times  near  the  axis,  the  radial  outflow  driven  by  Er  Is  turned  forward  by  the 
magnetic  field  and  the  net  plasma  current  density  Is  forward.  The  amount  of 
current  Involved,  however.  Is  Insignificant  In  magnitude  compared  to  the  beam 
current. 


When  the  high-energy  electron  group  from  the  3-group  model  Is  examined 
separately.  It  shows  quite  large  effects.  The  z-component  of  plasma  current 
for  the  high  group  Is  shown  at  various  times  In  Figure  2.15(a).  Beyond  about 
1.5  Bennett  radii  it  Is  always  directed  forward  with  the  beam.  Nearer  the  axis, 
the  strong  Ez  field  keeps  the  current  going  backwards,  although  by  10  ns  there 
is  relatively  little  backward-moving  current  In  the  high  group.  A  small  part 
of  this  forward-going  current  is  due  to  the  fact  that  high-group  electrons 
produced  directly  by  the  beam  are  Injected  with  forward  velocity.  However, 
most  of  It  is  due  to  the  Hall  force  acting  on  the  radial  out-flow. 

The  effect  of  Hall  current  In  th.e  high  energy  electron  group  Is  shown 

In  Figure  2.15(b).  At  1  ns,  the  forward-going  current  near  1  Bennett  radius 

2  2 

amount  to  almost  2  kA/cm,  compared  to  the  ^  5  kA/cm  of  low  group  plasma 
current  going  backwards  at  the  same  radial  position.  The  result  Is  a  very 
significant  reduction  of  the  total  return  current. 

It  Is  Important  to  note  that  at  late  times  £  >  0.5  ns,  the  radial  flow 
which  is  turned  forward  by  the  magnetic  field  Is  driven  not  primarily  by  Ep, 
but  by  the  radial  pressure  gradient  of  the  high-group  particles.  Yu  (Ref.  4) 
first  recognized  this  as  a  possible  Important  component  of  the  total  current. 
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Figure  2.15(a).  Current  Density  0  of  High  Energy 
Electrons  (p/pQ  *  .008). 
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Figure  2.15(b).  Effect  of  High  Energy  Electrons 
on  Total  Plasma  Current  Density  J 
at  1  ns  (p/p  *  .008). 
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From  the  equation  of  motion  It  Is  clear  that  a  density  gradient  Is  equivalent 
to  a  radial  electric  field  of  magnitude 

Er  *  -  T  -gvy—  volts/cm. 

For  T  'v  10  KeV  as  suggested  by  Yu  (and  used  for  the  high  group  In  these  cal¬ 
culations),  and  a  density  gradient  with  scale  <\<  1  cm,. .this  corresponds  to 
Er  'v*  10  KeV/cm.  Assuming  that  the  system  as  a  whole  approaches  a  steady  state 
solution  of  Maxwell's  equations,  must  -*•  0.  In  the  presence  of  the  Hall- 
force  terms,  the  radial  Er  adjusts  Itself  to  shut  off  the  net  current  driven 
by  v2B0,  Er,  and  the  pressure  gradient.  In  the  present  context  it  seems  that 
zero  net  radial  current  could  perhaps  be  achieved  by  a  balance  between  an  In¬ 
ward  flow  of  very  low  energy  electrons  and  an  outward  flow  of  a  much  smaller 
number  of  fast  electrons.  Because  of  the  differences  In  v  for  the  two  streams, 
the  corresponding  Hall  currents  would  not  cancel  exactly  In  this  case.  Under 
some  conditions  such  counterflows  may  be  limited  by  instabilities. 

2.5.5  Delta  Ray  Effects  at  Low  Density 

Although  the  steady  state  5 -ray  current  Is  Independent  of  air  density, 
the  time  required  for  the  current  to  reach  its  maximum  value  depends  strongly 
on  the  density.  There  is  also  a  significant  dependence  on  the  beam  energy  if 
only  relativistic  particles  (>  1  MeV)  are  included,  as  shown  in  Figure  2.16. 

The  pulse  length  used  In  this  calculation  was  100  ns.  (A  very  steep  drop  In 
6 -current  at  1  atmosphere  on  a  timescale  of  10  ns  beginning  just  before  the 
pulse  ends  Is  not  shown.)  It  Is  clear  that  at  densities  as  low  as  .01  atmosphere 
the  5-current  represents  a  very  small  Increment  to  the  beam  current  over  the 
first  10  ns  considered  In  the  calculations  presented  above.  Before  0.1  ns, 
the  5 -current  Is  comparable  to  the  plasma  current  In  the  z-directlon,  but 
both  are  very  much  smaller  than  the  beam  current.  It  seems  unlikely  that  the 
relativistic  part  of  the  beam  secondary  cascade  can  be  Imrt  *.'nt  at  early 
times  for  air  densities  as  low  as  .01  atomsphere.  However,  for  propagation 
In  a  density  channel,  the  interaction  with  the  highrdensity  walls  may  be  very 
Important. 
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Figure  2.16.  (Delta  Ray  Current)/ (Beam  Current)  at  Energy  >  1  MeV. 
Solid  Curves:  Ebeam  3  100  MeV 

Dashed  Curves:  Ebeam  3  10  MeV 
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2.5.6  Comparison  with  Ohm's  Law 

Both  the  1-  and  3-group  models  use  equations  of  motion  to  describe 
the  electron  flows.  In  general,  the  difference  between  the  currents  calculated 
this  way,  and  by  using  an  Ohm's  law  relationship  Involving  only  electron  density 
and  momentum  transfer  frequency,  Is  large  at  early  times.  The  duration  of  the 
deviation  from  Ohm's  law  decreases  as  the  gas  density  decreases.  Figures  2.17(a) 
and  (b)  show  the  actually-computed  plasma  current  compared  with  what  would  have 
been  obtained  from  an  ohmic  calculation  using  the  Instantaneous  values  of  Ne 
and  vm.  The  results  are  shown  for  the  1-group  calculation  at  p/pQ  *  .008,  but 
comparable  effects  are  present  for  the  3-group  case.  At  0.1  ns  the  ohmic  current 
shows  an  extremely  strong  Hall  effect  due  to  the  .large  radial  field;  the  current 
given  by  the  equation  of  motion  Is  very  much  smaller,  and  the  Hall  effect  does 
not  show  on  the  linear  scale.  At  1  ns,  the  two  methods  of  calculating  the 
current  give  virtually  identical  profiles,  as  shown  In  Figure  2.17(b). 

2.5.  7  Concluding  Remarks 

The  calculations  discussed  above  show  the  magnitude  of  the  effects  to 
be  expected  when  low-density  phenomena  are  taken  Into  account.  All  the  comnonly 
used  assumptions  listed  In  Section  2.2  above  have  significant  Impact  on  the 
results  at  p/pQ  *  .01.  (The  current  carried  by  relativistic  beam  secondaries 
Is  Important  at  low  densities  only  If  the  pulse  is  long  enough  for  the  full 
build-up  to  occur.) 
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Figure  2.17(a)  and  (b).  Comparison  of  Ohmic  Calculation  of 

Plasma  Current  Density  with  Equation 
of  Motion  Calculation. 


i 


f. 

-*  V.V; 


3.0  HERMES  HOSE  EXPERIMENTS 

During  the  period  01  September  1983  through  31  March  1983,  SAI  provided 
theoretical  support  to  Sandla  National  Laboratory,  Albuquerque  (SNLA)  to  deter¬ 
mine  the  relevance  of  the  HERMES  facility  to  address  key  uncertainties  In  high- 
current  beam  hose  stability.  The  theoretical  conclusion  was  quite  optimistic: 

With  proper  pulse  conditioning,  and  with  small  initial  hose  perturbations,  hose 
growth  should  be  small  enough  to  permit  adequate  calibration  of  existing  hose 
codes.  Unfortunately,  In  the  limited  amount  of  time,  spent  on  the  HERMES  machine. 

It  was  not  possible  to  solve  all  the  difficulties  associated  with  beam  Injection. 
Although  the  beam  radius  was  small  enough  to  be  of  Interest  (1-2  cm),  the  pulse 
was  often  Injected  off-axis  (probably  due  to  non-uniform  cathode  Ignition). 

Various  conditioning  cells  were  employed  with  some  success,  but  the  Initial  hose 
disturbance  was  always  too  large  to  compare  with  low-amplitude,  linearized  hose 
models.  Although  theory  predicts  that  hose  should  saturate  after  about  6  meters 
of  propagation  (with  HERMES  beam  parameters),  the  measured  hose  displacements 
after  only  50  cm  of  propagation  were  already  In  the  non-linear  regime. 

It  was  noted  by  SNLA  that  in  all  drift  chamber  pressures  a  "low- 
frequency"  hose  sweep  could  always  be  detected  (by  the  Ba  probes  and  the  Faraday 
cups)  with  frequencies  between  30  and  70  MHz.  This  is  in  qualitative  agreement 
with  the  high-frequency  cutroff  prediction  of  all  multi -component  hose  models. 

At  a  fixed  position  z  from  the  injection  point,  the  dominant  sweep  frequency 
as  the  pulse  passes  should  vary  inversely  with  the  magnitude  of  the  instan¬ 
taneous,  magnetic-dipole  relaxation  time  The  variation  of  with  time,  chamber 
pressure,  and  the  distance  z  as  predicted  by  the  PHLAP  code  is  in  qualitative 
agreement  with  the  observed  variations  In  hose-sweep  frequencies. 

The  SAI  code  PHLAP  (Ref.  9)  was  employed  to  predict  the  hose  growth 
of  the  HERMES  pulse  for  various  pressures  and  initial  conditions.  The  beam 
current  profile  was  taken  directly  from  measured  values  and  Is  plotted  In  Figure 
3.1  as  a  function  of  time  c  (nsl  at  the  fixed  position  z  *  0.  The  retarded  time 
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variable  €  measures  the  distance  back  from  the  pulse  head.  Although  HERMES 
produces  pulses  with  peak  current  greater  than  100  kA,  pulse  conditioning  reduces 
this  to  40  kA.  The  beam  energy  is  assumed  to  vary  linearly  from  3  MeV  at  5  *  0 
to  10  MeV  at  %  -  40  ns  and  is  held  at  10  MeV  for  5  >  40  ns.  The  initial  emittance 
is  chosen  to  match  the  pulse  at  5  *  35  ns  In  air  at  200  Torr  with  a  Bennett  radius 
of  3  cm. 

In  Figure  3.2,  the  displacement  of  the  center  of  charge,  60  ns  behind 
the  head,  is  plotted  as  a  function  of  propagation  distance  z.  The  initial  dipole 
perturbation  was  random  noise  with  a  maximum  amplitude  of  unity.  (The  PHLAP 
dipole  model  is  completely  linear.)  The  hose  growth  saturates  at  about  z  =  6  m 
with  a  maximum  gain  of  46  and  then  commences  to  damp-out. 

In  Figure  3.3,  the  maximum  saturated  gain  is  plotted  for  all  positions 
behind  the  pulse  head  at  two  different  propagation  chamber  pressures;  50  and  200 
Torr.  The  saturation  of  maximum  hose  growth  in  the  pulse  tail  is  an  important 
high-current  prediction  and  hence,  could  in  principle  be  tested  with  HERMES  or 
other  existing  diode  machines.  In  addition,  the  scaling  of  hose  growth  with 
pressure  could  also  be  calibrated  in  near  term  experiments. 
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